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Abstract 



We show the evolution of a gluon plasma towards equilibrium starting at some early 
moment when the momentum distribution in the central region is moment aneously 
isotropic. Using HIJING results for Au+Au collision as initial input, we consider ther- 
malization and chemical equilibration simultaneously at both LHC and RHIC energies. 
Thermalization is shown to be driven chiefly by inelastic process in our scenario contra- 
dicting common assumption that this is the role of elastic process. We argue that only 
the inelastic dominancy depends on the initial conditions but not the dominance itself. 
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1 Introduction 



In future heavy ion collision experiments at the Relativistic Heavy Ion Collider (RHIC) 
at Brookhaven and at the Large Hadron Collider (LHC) at CERN, one hopes to create 
a new state of matter, the quark-gluon plasma. In order to clearly distinguish this from 
the possible alternative of a hadron gas, one would like to know the state of the plasma 
in its different stages throughout its lifetime. With this knowledge, the conditions under 
which different possible signatures of the plasma are produced and the time at which 
they are produced can be compared with those of the hadron gas |3J. It is precisely the 
difference of these conditions and the time dependence which could permit a possible 
distinction from a hadron gas and a clear identification of a quark-gluon plasma. There 
have already been numerous works in the direction of the possible signatures and also 
a number of investigations on the various time scales and stages of the plasma. In this 
work, we will focus on the latter. 

In the following, we will not attempt to match the comprehensiveness and details of 
the semi-classical parton cascade model P, [| which models the evolution in relativistic 
heavy ion collisions by tracing the phase space history of the particles and combining 
perturbative QCD with relativistic transport theory. Our approach will base on rather 
simple assumptions and initial conditions which have already been made used of in the 
previous studies of thermalization and parton chemical equilibration @, 0, 

In many previous studies |], [|, ||, |[ [| of relativistic heavy ion collisions, a con- 
venient assumption and starting point is kinetic equilibration is very rapid < 1 fm and 
hydrodynamic expansion is well underway. As has been pointed out in [[j], [HJ , collisions 
are not sufficiently rapid at early times to maintain thermal equilibrium in an expanding 
plasma, so the system is more likely to be in a situation between free streaming, where 
there is no collision or the collision effects are completely negligible, and hydrodynamic 
expansion when local kinetic equilibrium has been achieved. Various studies have been 
made in determining the relation between collision time 9 and the degree of equilibration 
of the plasma Jf], [TTj, based on relaxation time approximation and assumed simple 
power behaviour of 9 ~ r p . In this work, we will also look at the collision time and 
thermalization but we do not assume a priori any form for 9, rather it will be determined 
entirely from the interactions. Our starting point will be different from many previous 
works in that we start somewhere between free streaming and hydrodynamic expansion. 
Some initial steps in this direction have been taken by Heiselberg and Wang in J7| but 
they only considered elastic interactions and only in the time range tq < r <C 9, where 
To is the initial time which will be further explained below. Following their initial line of 
approach, we attempt to give a more complete treatment. We will do this under certain 
assumptions and approximations which will be explained in the subsequent sections. 
The method that will be shown here allows the treatment of the out of equilibrium 
plasma with comparative ease. We will illustrate this approach by applying it to the 
simpler case of a pure gluon plasma. Thermal as well as chemical equilibration of the 
gluon plasma will be investigated. That of a quark-gluon plasma will be dealt with in a 
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subsequent paper [JL3 



For the initial conditions, we will use results obtained from the Monte Carlo HIJING 
model [14, 15, 16| as have been used in 0] and later improved in || [(J but we will interpret 
them in a slightly different way so that our initial temperatures T and fugacities Iq are 
not the same as theirs although we use the same initial gluon energy densities eo and 
number densities no- 

In principle, one could vary the initial conditions and look at the resulting evolution 
of the plasma as have been done in || since HIJING predictions for the energy densities 
and number densities are in the low end in comparison with, for example, those from 
the parton cascade model |§ but we will not do that here. We would rather concentrate 
on one set of specific initial conditions each at RHIC and LHC energies and look at the 
thermalization, isotropy of momentum distribution and the time dependence of fugacity, 
collision time 9 and a modified form of the non-equilibrium colour screening length. We 
show, in particular, inelastic process dominates over the elastic one in the approach 
towards equilibrium and this is independent of initial conditions. The degree of this 
dominance, however, is dependent on initial conditions. 

The paper is organized as follows. In Sect. |2], the various assumptions, approxima- 
tions and approach will be described. Sect. is devoted to calculating rates for the 
interactions that we will consider and the set of equations to use for chemical equili- 
bration will be given. Initial conditions used and the method of computation will be 
outlined in Sect. |. Lastly, the results of the evolution of the gluon plasma will be 
shown and discussed in Sect. ||. 



2 Boltzmann Equation and Relaxation Time Ap- 
proximation 

In a many-body system, it is expedient to use particle distributions /(p,r,r) as well 
as collective variables energy density e, number density n, entropy density s etc. or 
at least local collective variables to describe the spatial as well as temporal variation 
of the system. This is what we will do here but we concentrate mainly in the central 
collision region which we assume to be essentially uniform in space. Ideally, one would 
like to describe f(p,r) by a set of fully relativistic and quantum transport equations 
derived from first principle of QCD [T7], [Tj| since a quark-gluon plasma is inherently 
a quantum system. But unfortunately, this has not yet been realized although various 
advances have been made in that direction based on constructing Wigner operators 
from the fields and deriving their equations of motion from the field equations JT^, [2(J . 
The difficulty is partly due to the complexity of the equations of motion for the Wigner 
operators which requires, in the end, one to make a semi-classical expansion in powers of 
h and partly due to the gauge degrees of freedom which make it hard to identify a gluon 



distribution operator [ZIJ. These gauge degrees of freedom of QCD have proved to be 
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a stumbling block in identifying a Wigner operator as the gluon distribution operators. 
Recent works [^2], [23| suggest that this latter problem could be overcome by choosing 
the radial gauge but one has yet to derive the collision terms. 

Under these circumstances, one has to make do with relativistic but semi-classical 
approach such as using Boltzmann equation 



d 



d 



/(p,r,t) = C(p,r,t), 



^dt dr 

and combining it with collision term C obtained from QCD. Assuming one-dimensional 
expansion along the z-axis |24| and boost invariance or rapidity independence of the 
central region, Baym [11] showed that Eq. ([!]) can be rewritten as 

( ~ =C(P±,Pz,T). 

IT| , we use the relaxation time approximation for the collision term 

f(p,-r) ~ feq(P,T) 



As in 



C(P,- 



(2) 



6(r) 



which enables us to write down a solution to the Boltzmann equation Eq. 

/ = fo(p±,PzT/r )e- x + 



dx'e x x f t 



eq v 



pl + ( Pz r/r'y,T eq (r')) 



where 



x(t) = / dr'/e{r') 



(3) 



(4) 



(5) 



given some initial distribution / at r . fo(p, r) is the free streaming, C — 0, solution 
to Eq. (0). Our solution differs from that in 0. In their case, they considered only 
elastic interactions and so had to keep a chemical potential /i(r') in f eq . We will not 
restrict ourselves to only these interactions. On the contrary, we will use the interaction 
gg < > ggg to determine 9{r). T eq (r) in f eq is determined from energy conservation 

e(r) = e eq (r) = a 2 T* q (r) , (6) 

with ci2 = 87r 2 /15. It is the momentaneous target temperature that the gluon plasma is 



striving to reach. The associated energy conservation equation is |TT 

de 



+ 



e +Pl 







dr r 

and pl is the longitudinal pressure in the z-direction given by 



(7) 



Pl{t) 



V 



P Pz 



f{p±,Pz,r) 



(2vr) 3 p 

with v = 2 x 8 = 16 for the multiplicity of gluons. Along the same line, we can also 
defined px, the transverse pressure by replacing p z by p x or p y in Eq. (|8J). These pressures 
will be used to check the degree of anisotropy of momentum distribution in the central 
region of the plasma. 
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3 Rate Equations 



a 



To consider chemical equilibration, we require the rate equations for gluon creation or 
destruction. From [25| and Eq. (Q), we have the following rate equation for the gluon 
number density 

and using Eq. (0) and Eq. ([!]), we also have 

m v _ M 7 ") - n(r) 

H(r) ^ , (10) 

where n eq {r) = aiT^ q {r) and a x = 16£(3)/7r 2 . The a in Eq. (|j) is a numerical factor to 
take into account of multiplicity and symmetry of the permutations of the gluons in C . 
The two equations Eq. (Q) and Eq. ( |T0"D allow us to determined 6{r) at any instance in 
time. 

For the r.h.s. of Eq. (|9]), we take this to be the simplest gluon multiplication process 
99 i * 999 

x[/i/ 2 (i + / 3 )(i + + h) - 73/4/5(1 + A)(i + / 2 )] 

x9(A-r QCD ). (11) 

The various parts of this formula need some explanation, we do this one by one. Here 
we only consider small angle scattering and without loss of generality p§ is taken to be 
a soft gluon radiation around zero rapidity in the centre of momentum frame so that we 
can use Bertsch and Gunion formula [^6[ for the amplitude. In that frame of the two 
incident gluons, where they approach each other on the z-axis, this is 

\ M I 2 = S - ) ( ) (12) 

gg ^ 9 " 1 W 2 -l(q 2 L + m 2 ) ) 2 A k ^( k± _ q± )2 + m 2 ) ]; K^) 

where and q_i_ are the perpendicular component of ps and that of the momentum 
transfer in that frame respectively. iV = 3 is for SU(3) of QCD and s is the total 
centre of mass energy. The first part on the r.h.s. within round brackets is the squared 
modulus of the small angle elastic scattering amplitude for gg — > gg and the second part 
is essentially the soft gluon emission spectrum. Since the centre of momentum frame 
of any two incident gluons are very unlikely to be the same as the plasma rest frame, 
we need a general formula for this amplitude. In the plasma rest frame, one simply 
substitues 

(k ± -ci ± ) 2 = Ap r p 3 p2-p3/s (13) 
qi = 4pi-p 4 p 2 -p 4 /s (14) 
k 2 L = Ap 1 -p 5 p 2 -p 5 /s (15) 
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in Eq. (U). 

We have infrared regularized the squared modulus of the amplitude Eq. ([12]) with 
the Debye screening mass m 2 D which is time dependent in the present case. We take it 
to be defined at leading order in a s by 



Note that the screening mass should be directional dependent in general |2q] but 
rather than keeping tracks of the directions of the gluons, we simplified this by removing 
the directional dependence. The present form can be thought of as an averaged screening 
mass between the transverse and longitudinal screening. We will explain this in some 
details elsewhere [13|. Here / is not of the Bose-Einstein form when the plasma is out of 



equilibrium but is given by Eq. (f|). It reduces, of course, to the familiar m 2 D = Ana s T 2 
when in equilibrium. This provides us with a general Debye screening mass to screen 
off infrared divergences even in out of equilibrium situation which is the subject of this 
paper. 

As in [El I5L H|, we use the approach of Gyulassy and Wang |2§|, [5(J to deal with 
the Landau-Pomeranchuk-Migdal (LPM) suppression of soft gluon radiations due to 
multiple scattering of the parent gluon in a medium. This is taken into account in the 



form of the ^-function in Eq. (|TT|). The basis of this approach is to include only the 
non-suppressed Bethe-Heitler limit but the suppressed emission between that and the 
factorization limit is completely dropped. Ideally, one would also like to include these 
emissions. The ^-function requires that the mean free path, A, of the gluon which we 
take from the small angle gg elastic scattering cross-section j|, ]j| , 



dair 99 9 2va\ 



dqi 4 (qi + m 2 D ) 2 

. , p/4 , 9 da 9 p r gg 9na 2 sn 
A = n dcf A ' 



(17) 



/o ""*" L dq]_ Sm 2 D (s / A + m 2 D ) 

to be larger than the formation time tqcd — ^o/k^ of the soft gluon. This LPM effect 
on the radiative energy loss of a fast energetic parton previously derived in has 



been improved and rederived recently in [31] by including interference graphs where the 



emitted gluon also undergoes multiple scattering with the scattering centres. However, 
the derived soft radiation LPM spectrum has no dependence on the direction of the 
emitted gluon. This has to be integrated out for the radiation density to make sense. 
As such, their result cannot be used in the present problem. One will have to be 
contented to include only the Bethe-Heitler non-suppressed emission therefore we use 
the same 6 function as in [Ej ^ ||. 

At this point, a remark is in order concerning large angle scattering non-soft gluon 
emission. In H, Shuryak and Xiong found that large angle scattering is also important 
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and higher order tree diagrams of maximal helicity violating amplitudes, the so called 
Kunszt-Stirling corrected Park- Taylor formula |33| for gluon multiplication are not 
negligible. In performing their calculation, LPM suppression has been totally neglected. 
As far as the author is aware, possible interference effect on non-soft gluon emission 
in multiple scattering has not been worked out. May be this is negligible for non-soft 
radiation if one is permitted to generalize the language of the soft emission to the non- 
soft case. That is by arguing that the gluon mean free path is approximately the same 
or larger than the coherence length of these radiations, but the latter is not defined 
for non-soft radiation. Facing the choice of giving up LPM entirely or including large 
angle scatterings, we choose the first alternative. The latter will be investigated in a 
future work where some form of LPM or more general interference effect due to multiple 
scattering will have to be included. 



4 Initial Conditions for the Computations 

As in ^j, H [7|, we consider the following scenario. Initially, the two incoming heavy 
ions which have been accelerated to very high energies, 200 GeV/nucleon at RHIC and 
up to 6.3 TeV/nucleon at LHC, along the z-axis collide, penetrate each other leaving 
behind in the central region a gas of secondary partons produced via minijets [34 . 
Since most energies are in the longitudinal direction, the momentum distribution is 
mainly along the beam direction. As mentioned in the introduction, collision time is 
larger than the expansion time so initially the partons are more likely to stream freely 
in the longitudinal direction gradually suppressing the momentum distribution in that 
direction in the central region. One can understand this by looking at the form of fo in 
Eq. (|2"2"D. At certain moment r , longitudinal momentum distribution will be reduced 
to such an extent that it becomes the same as the transverse momentum distribution, 
or in other words, the momentum distribution becomes isotropic. Assuming a thermal 
distribution at this time r fixes fo(p,r Q ). The initial temperature T and fugacity l 
can be determined from the gluon energy density and number density. These values 
will be taken from the improved version || of Q but we obtain the temperature and 
fugacity in a slightly different way. In |4], ||, they used the fugacity factorized thermal 
equilibrium form of /, i.e. 

/(p,Q = n im i (19) 
exp(|p|/T) - 1 

for all values of I. As was pointed out in their work, this is a good approximation only 
when I « 1. Here we keep the original form of / at least whenever it is in thermal 
equilibrium. We get the initial temperature T and fugacity lo at r by observing that 
Iq is small initially, at least in HIJING results, when we can approximate / by the 
Boltzmann form / = I e~ p / T so the energy density €q and number density tiq are given 
by 

e = 3 v l T 4 /tt 2 and n = ul T 3 /vr 2 , (20) 
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respectively. The extracted To and Iq for Au+Au collision at LHC and at RHIC are 
shown in Table 1. Knowing 

/o(P,T ) = 7j k ' r , (21) 

exp(|p|/T ) - Z 

we can deduce the time dependent form for it, since it is the free streaming solution to 
Eq. §10 

fo(p±,Pz,r) = l ° . (22) 

exp( v /pi+^ r 2 /r 2 /To) _ /q 

Using this, we can calculate the rate gg < — > ggg at tq and hence 9q. With these initial 
values, we start the evolution of the plasma. 

The evolution of the gluon plasma is performed as follows. Phase space integrals 
for calculating the gluon multiplication rate are performed by an iterative and adaptive 
Monte Carlo method VEGAS [|3~5]| . All time integrals are calculated by discretizing time 



r so that any function G(t) depending on a time integral of another function g(r, r') is 
re-expressed as 

/ T n— 1 

dr'g(r, r') => G{r n ) ^^Ar g(r n , r t ) , (23) 

D i=0 

where Tj = r + iAr and Ar is some small chosen time step. Then G(r n ) is determined 
entirely by variables of the preceding time steps from r to r n -\. Given all the variables 
at to, the same set of variables at T\ can be determined. This procedure is repeated for 
obtaining the variables at T2 from those at To and T\ and so on. The set of variables to 
determine are e(r) for obtaining T eq (r) which enables the determination of m 2 D (r). This 
modified screening mass is required for calculating R(r) which together with n(r) and 
T eq (r) allow us to obtain 6{t). Both 6{r) and T eg (r) are essential for finding /(p,r) 
which is itself needed in R(t) and for the calculations of the various collective variables. 



5 Equilibration of the Gluon Plasma 

The results of the evolution towards equilibrium of a gluon plasma for two sets of initial 
conditions for Au+Au collisions at LHC and at RHIC are presented in Fig. 1, 2 (LHC) 
and Fig. 3, 4 (RHIC). These are results calculated using a s = 0.3 and Ar = 0.05 
fm/c. Using smaller values of Ar, for example, Ar = 0.02 fm/c shifts the various 
values by amounts ranging up to 6.25%. Since we are not going after high numerical 
accuracy, we will settle for the present value. In any case, there is no qualitative change 
by varying Ar. Unless explicitly stated otherwise, the dashed curves in most figures 
are the equilibrium values associated with f eq in Eq. (|3|). The solid curves are the 
actual results. One sees that the two curves approach each other with increasing time. 
How close the solid and dashed curves are is a measure of how close the plasma is to 
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equilibrium. The computation is stopped when the estimated temperature falls below 
200 MeV. This and the fugacity are estimated by using the expression for the near full 
equilibrium energy and number density 

e = a 2 lT 4 and n = a x I T 3 (24) 

which are valid when I is close to one, so that 

T es t = ai e/ {a 2 n) and l est = a\ n 4 / (a\ e 3 ) . (25) 

Since we have used energy conservation e(r) = e eq (r) at all time so in Fig. 1 and 3 
(a) where we show the evolution of the energy densities, the solid and dashed curves 
coincide. Comparing gluon densities in Fig. 1 and 3 (b), the estimated temperatures 
in Fig. 1 and 3 (c), the entropy densities in Fig. 1 and 3 (e) and the squared of the 
modified screening mass in Fig. 2 and 4 (b), it is apparent that at LHC, the solid and 
dashed curves are much closer to each other at the end of the calculation than they are 
at RHIC. Those at LHC are actually lying on top of each other or at least very close. 
The estimated fugacity l est at the end of the calculation is about 0.94 at LHC and 0.40 
at RHIC. The time evolution of these and the estimated temperatures are shown in 
Fig. 1 and 3 (c). The dot-dashed curves are the fugacity evolutions. Note that these 
estimated values can only be viewed as indicators. It will be seen that the plasma have 
not quite completely thermalized at the end of the calculations. 

The entropy density evolutions s(r) shown in Fig. 1 (e) for LHC and in Fig. 3 (e) 
for RHIC are calculated from PS| 

s(r) = -vj ^s{/(P, t) In/ (p, r) - (1 + /(p, r)) ln(l + /(p, r))} . (26) 

The dashed lines are s = Ae/3T eq , the ideal gas form of the equilibrium entropy. To 
look at the actual time dependence of the generated entropy, we plotted the product of 
the entropy density and r, st, which is proportional to the total entropy in the central 
region in Fig. 1 and 3 (f). At LHC, the entropy no longer increases or stopped being 
produced at the end at r = 8.75 fm/c, whereas at RHIC, it has not yet stopped to 
increase at r = 3.55 fm/c when the estimated temperature falls below 200 MeV. 

To check for isotropy in momentum distribution, we plotted the ratio of longitudinal 
to transverse pressure in Fig. 1 and 3 (d). These pressures are as defined in Sect. |J. 
The solid curve in each case is this ratio while the dashed curve is the ratio of a third 
of the energy density to the transverse pressure for comparison. The plasma starts off 
in an isotropic momentum configuration and becomes anisotropic due to the tendency 
of the plasma to free stream in the z-direction. This anisotropy attains its maximum 
at about 1.35 fm/c at LHC and 2.2 fm/c at RHIC. Since the plasma is approaching 
equilibrium, both Pl/pt and e/3pT must approach 1.0 with increasing time. The figures 
show the plasma has not quite regained isotropy in momentum distribution although it 
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is much closer at LHC than at RHIC. Note that pi < e/3 except at the start therefore 
less work is required for the expansion than in the case of hydrodynamic expansion 
when p L = e/3. The consequence is the energy density drops less fast with increasing 
t. Note that this does not necessarily mean that the plasma can last longer than if it 
underwent hydrodynamic expansion. Particle production rate must also be compared 
because although this does not change the energy density, it does reduce the temperature 
whenever the concept of a temperature makes sense. If the rate is not so rapid that it 
over compensates for the reduction in pl, then indeed the plasma can last longer than 
that undergoing hydrodynamic expansion before phase transition takes place. 

To look at the relative importance of free streaming and collisions, we plotted in 
Fig. 2 and 4 (a), the modified screening mass squared (solid curves) and that of free 
streaming when collisions have been switched off. The modified screening mass squared 
is 



m 2 D (r) 



^'(ro/r) ml e "*M + Ana s e~^ jf ^X^'M T^r') e°W (27) 



where 



sin 1 a/1 — r 2 



Kr) = rr~ o (28) 



and that of free streaming is obtained by letting 9 — > oo. The latter is a decreasing 
function of time so free streaming increases the screening length inside the plasma since 
the gluon density decreases with r -1 . Whereas collisions tend to shorten the screening 
length or increase the screening mass due to the equilibrium value is several times higher. 
We see that even at the beginning, collision effect is already very strong especially at 
LHC where 9o ~ 0.84 when t ~ 0.5 so the plasma starts very much in the mixed phase 
of free streaming and hydrodynamic expansion. This is less so at RHIC with 6$ ~ 2.7 
when ro ~ 0.7 but the collision effect is still strong. In both cases, the screening mass 
increases first with time (more rapid increase at LHC than at RHIC) before it decreases 
eventually because the plasma is undergoing expansion and therefore cooling. Fig. 2 
and 4 (b) show the convergence of m 2 D (r) with the equilibrium value with increasing r. 
Again, we see a better convergence at LHC than at RHIC. 

Another evidence for the strong collision effects is the time dependence of the collision 
time 6{t) plotted in Fig. 2 and 4 (c) (These curves are relatively rough but because 
of the time required for the computations, we do not attempt to get smoother curves. 
In any case, how 9 is evolving with time is clear.). One expects free streaming to drive 
the system further out of equilibrium and hence increases 9. Fig. (c)'s show quite the 
opposite behaviour of a clear initial drop in 9 in each case. A much steeper drop at 
RHIC than at LHC is seen. These initial drops are followed by relatively slow increase 
with respect to r. These behaviours are essential for the plasma to be able to reach 
equilibrium i.e. they allow r to catch up with 9 and overtake it. This overtaking of 9 
by r is much more dramatic at RHIC since 9 is about 4 times larger than r whereas 
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at LHC, 9q/tq is only about 1.6. The final slow increase in 9 is expected from previous 



calculations of relaxation time |38|. They found in general when the system is 
near equilibrium a dependence r re i ax ~ 1/T on the temperature which means a time 
dependence of r re i ax ~ r 1//3 . In our case, we can parametrize the rate R(t) near full 
equilibrium by 

R(t) « (I-/) cT 4 = (1-/) c'aiT 4 (29) 
where c and d are some constants. Using Eq. (f24"|), we get 

An = n eg -n = a x T 3 (Z 3/4 - I) . (30) 

From the above two equations, one obtains the near equilibrium collision time #z_>i(t) 

/3/4 

ft-iOO = cT( i + / i/ 2 ) ( i + / i/4) • ( 31 ) 

We have plotted this in Fig. 2 (c) (dashed curve) using T est and l est but not in Fig. 
4 (c) since at RHIC, the plasma is still not yet near equilibrium. The constant d 
determined from the rate is roughly d ~ 0.16377. One sees that the solid and dashed 
curve converges at large time. One could also fit a r 1//3 curve on the same plot (not 
shown for clarity) which would show 9 does have the right behaviour for a equilibrated 
or nearly equilibrated plasma. 

Although our 9 has a more complicated time dependence than simple power be- 
haviour which was heuristically assumed in previous works [0, [11], [L2| to enable simple 



analysis, it does satisfy the generalized condition of Heiselberg and Wang for equili- 
bration at all time. Their analysis for 9 ~ r p showed thermalization would be achieved 
for p < 1 otherwise 9 would increase faster than r or collision would always be slower 
than the expansion. For general 9, this can be reworded in terms of x(r) defined in 
Eq. (]5D as x — > oo as r — > oo for thermalization, if x — > x max < oo as r — > oo, then 
thermalization cannot be realized. 

Finally, a word must be said about elastic collisions gg — > gg which have almost been 
left out entirely since they drop out of Eq. (|S|) (their only role so far is in determining 
the cross-section). In 0, by considering only elastic processes, 9 has been determined 
from the collision entropy density rate 



ds 



coll 



-"J ' Wf C{ *' T) ^TTf) (32) 



where a' is again some numerical factor for multiplicity etc.. They found that 9 had 
only logarithmic time dependence. Their result, which is approximately valid in the 
range 2r £ t 9 and after being readjusted to our initial conditions, is 

^ 1.55^/0^0 In ( 3?rr )ln(— ) . (33) 
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Using t between 2tq £ r <, 9/2, we find 9 varies between 4.28-3.53 fm/c from r ~ 1.40 
fm/c to r ~ 1.76 fm/c at RHIC, which is about 2.0 times higher than our 9 in the 
same time range. At LHC, the above time range does not exist and so Eq. (33) cannot 
be applied. This shows that inelastic process gg < — > ggg seems to be more efficient 
for thermalization than elastic gg — > gg at least at RHIC. To verify this, we have 
checked using Eq. ( p2[ ) how important is the entropy generated by elastic scattering 
when compared to that by inelastic scattering both at RHIC and at LHC. The result is 
small. We have to conclude that in our scenario the latter is the chief driving force for 
thermalization which contradicts the common assumption that thermalization is mainly 
the role of elastic processes. This conclusion is perhaps not too surprising since we have 
started in a momentaneously thermalized plasma with some small initial fugacities when 
all elastic processes are "out of action" . Our figures show that the plasma quickly reaches 
the point of greatest departure from equilibrium before heading back progressively. In 
order for elastic scatterings to play a more prominent role, there need to be some means 
for an explosive rise in rate before that pointQ has been reached. This means is, however, 
missing. It must be said that although our conclusion on the dominance of inelastic over 
elastic process is independent of the initial conditions, the degree of dominance is not. 
Had we started with some much larger initial fugacities, say lo ~ 1, then the initial 
gluon multiplication rate R(r = r ) would be greatly reduced and this will affect the 
dominancy. This can be understood in the following way. Since one cannot perturb an 
equilibrated gluon plasma kinetically without doing so chemically at the same time [J any 
perturbation will always drive a plasma in equilibrium further off chemical than kinetic 
equilibrium. In our case, the perturbation is the one-dimensional expansion and it is 
not obvious that the plasma has been chemically perturbed unlike the perturbation in 
the kinetic direction. This is so because the chemical perturbation due to the expansion 
is too small when compared to that with which we imposed on the system by having 
small Zq's. For Iq ~ 1, the rate in Fig. 2 and 4 (d) will not be monotonically decreasing 
with increasing r. They will instead first rise to a peak before coming back down again 
towards zero. A typical sign that the system is going out of equilibrium before making 



the return QT3J . In this case, clear chemical along with kinetic equilibrium are further 
upset by the expansion. Numerically, this has been verified with Iq ~ 0.8. In this case, 
the relative importance of the two types of process is reduced but inelastic remains the 
more important one. To summarize, it is because of the fact that chemical equilibrium is 
always further away than kinetic equilibrium, inelastic process will always be significant 
in relation to elastic process in equilibration. A small Iq will give it a bigger advantage 
or increased importance than a lo ~ 1 



2 We can reasonably expect the elastic scattering rate maximizes at this point. This is within the 
time range in which we compared #'s above. 

3 The reverse is, however, not true. One can see this by finding a small perturbation that will hold 
chemical but not kinetic equilibrium. 
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Figure Captions 



Fig. 1 Results of the evolution of a gluon plasma towards equilibrium at LHC energies. 
The figures are (a) energy density e, (b) number density n (solid curve), (c) esti- 
mated temperature T est (solid curve) and fugacity l est (dot-dashed curve), (d) the 
ratio of longitudinal to transverse pressure Pl/pt (solid curve) and that of a third 
of the energy density to the transverse pressure e/3pT (dashed curve), (e) entropy 
density s, (f) the product of entropy density and time sr. The dashed curves in 
(b), in (c) and in (e) are the equilibrium number density, temperature and the 
ideal gas entropy s = 4e/3T eq associated with f eg respectively. 

Fig. 2 Results of the evolution of a gluon plasma towards equilibrium at LHC energies. 
The figures are (a) the modified screening mass squared m 2 D (solid curve) compared 
with the free streaming case (dashed curve), (b) m 2 D (solid curve) together with 
that in equilibrium (dashed curve), (c) collision time 9 and (d) the rate R for 

gg < > ggg. In this case only (not in Fig. 4), 9i^i is also shown in (c) (dashed 

curve) . 

Fig. 3 Results of the evolution of a gluon plasma towards equilibrium at RHIC energies. 
The figures and curves are organized in the same way as those in Fig. 1. 

Fig. 4 Results of the evolution of a gluon plasma towards equilibrium at RHIC energies. 
The figures and curves are organized in the same way as those in Fig. 2. 



Table 1 



Initial Conditions 




RHIC 


LHC 


r (fm/c) 
e (GeV/fm 3 ) 

n (far 3 ) 

T (GeV) 
k 

0o (fm/c) 


0.7 
3.2 
2.15 
0.50 
0.08 
2.69 


0.5 
40.0 
18.0 
0.74 
0.21 
0.84 



TABLE 1. Initial conditions for the evolution at RHIC and at LHC 
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FIG. 1, LHC 






FIG. 2, LHC 




FIG. 3, RHIC 





FIG. 4, RHIC 
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